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The  isolated  effect  of  02(a* 1Ag)  on  the  propagation  of  C2H4  lifted  flames  was  studied  at  reduced  pressures 
(3.61  kPa  and  6.73  kPa).  The  02(a1Ag)  was  produced  in  a  microwave  discharge  plasma  and  was  isolated 
from  0  and  03  by  NO  addition  to  the  plasma  afterglow  in  a  flow  residence  time  on  the  order  of  1  s.  The 
concentrations  of  02(a1Ag)  and  03  were  measured  quantitatively  through  absorption  by  sensitive  off-axis 
integrated-cavity-output  spectroscopy  and  one-pass  line-of-sight  absorption,  respectively.  Under  these 
conditions,  it  was  found  that  02(a1Ag)  enhanced  the  propagation  speed  of  C2H4  lifted  flames.  Comparison 
with  the  results  of  enhancement  by  03  found  in  part  I  of  this  investigation  provided  an  estimation  of  2-3% 
of  flame  speed  enhancement  for  5500  ppm  of  02(a1Ag)  addition  from  the  plasma.  Numerical  simulation 
results  using  the  current  kinetic  model  of  02(a1Ag)  over-predicts  the  flame  propagation  enhancement 
found  in  the  experiments.  However,  the  inclusion  of  collisional  quenching  rate  estimations  of  02(a1Ag) 
by  C2H4  mitigated  the  over-prediction.  The  present  isolated  experimental  results  of  the  enhancement 
of  a  hydrocarbon  fueled  flame  by  02(a1Ag),  along  with  kinetic  modeling  results  suggest  that  further  stud¬ 
ies  of  CnHm  +  02(a1Ag)  collisional  and  reactive  quenching  are  required  in  order  to  correctly  predict  com¬ 
bustion  enhancement  by  02(a1Ag).  The  present  experimental  results  will  have  a  direct  impact  on  the 
development  of  elementary  reaction  rates  with  02(a1Ag)  at  flame  conditions  to  establish  detailed 
plasma-flame  kinetic  mechanisms. 

©  2010  The  Combustion  Institute.  Published  by  Elsevier  Inc.  All  rights  reserved. 


1.  Introduction 

The  development  of  residence  time  restricted  systems,  such  as 
high-speed  air-breathing  propulsion  vehicles,  has  created  the  chal¬ 
lenges  of  achieving  efficient  and  reliable  ignition,  flame  propaga¬ 
tion,  and  flame  stabilization.  Supersonic  combustion,  for 
example,  constrains  the  flow  residence  time  to  millisecond  time 
scales  in  the  combustor,  therefore  limiting  the  time  for  chemical 
reaction.  In  this  case  the  flow  residence  time  is  comparable  to 
the  chemical  reaction  time  and  becomes  a  more  significant  prob¬ 
lem  when  using  hydrocarbon  based  fuels  instead  of  hydrogen  be¬ 
cause  of  slower  initiation  times.  Therefore,  the  development  of 
new  methods  to  decrease  ignition  delay  times  and  increase  flame 
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stabilization  and  flame  propagation  rates  are  critical  for  the  devel¬ 
opment  of  vehicles  capable  of  hypersonic  flight. 

There  are  two  methodologies  to  overcome  the  characteristic 
time  mismatch,  namely  increasing  the  flow  residence  time  or 
decreasing  the  chemical  reaction  time.  Modifying  the  flow  times 
through  geometric  changes  are  effective,  but  can  create  significant 
thermal  management  difficulties  and  pressure  losses.  Therefore, 
applying  new  techniques  to  increase  the  rate  of  fuel  oxidation 
and  overall  chemical  kinetics  has  attracted  much  attention. 

One  of  the  potential  solutions  to  increase  the  rate  of  combustion 
and  chemical  reactions  is  the  application  of  plasma  activation. 
Plasma-assisted  combustion  produces  elevated  temperatures,  rad¬ 
icals,  excited  species,  ions,  and  electrons  that  have  the  possibility 
to  increase  the  rate  of  fuel  oxidation.  Because  of  the  significant 
promise  of  plasma,  much  research  has  been  performed  using  a 
variety  of  plasma  discharge  systems  including  plasma  torches/jets 
[1-4],  gliding  arc  discharges  [5-7],  fast  ionization  waves  [8,9],  and 
nanosecond  repetitively  pulsed  discharges  [10,11],  as  well  as 
through  electric  field  interactions  [12-14],  microwave  discharges 
[15-17],  and  many  others.  The  investigations  have  shown  defini¬ 
tively  that  plasma  can  enhance  combustion  processes  with  de¬ 
creased  ignition  times  and  lower  ignition  temperatures 
[6,7,10,11,19,20],  increased  flame  propagation  [13,14,16-18],  and 
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Nomenclature 

A 

pre-exponential  factor 

Nozone 

number  density  of  ozone 

Ea 

activation  energy 

P 

pressure 

eV 

electron  volts 

ppm 

parts  per  million 

/ 

intensity  of  light  with  the  presence  of  ozone 

T 

temperature 

Io 

intensity  of  light  without  the  presence  of  ozone 

k 

reaction  rate  coefficient 

L 

absorption  path  length 

^ozone 

absorption  cross-section  of  ozone 

enhanced  flame  stabilization  [2,3,7,15].  It  has  been  stated  in  a  re¬ 
cent  review  of  plasma-assisted  combustion  that  “artificial  initia¬ 
tion  of  chemical  chains  by  low  temperature  plasma  of  gas 
discharges  does  exist,”  but  there  is  a  lack  of  data  from  controlled 
and  well-defined  experiments  [21].  Experiments  that  explore  the 
enhancement  of  combustion  processes  with  plasma  are  plagued 
with  highly  coupled  thermal  and  kinetic  interactions,  making  it 
difficult  to  isolate  the  key  species  and  the  various  enhancement 
mechanisms. 

Therefore,  one  of  the  major  challenges  associated  with  under¬ 
standing  the  plasma-flame  interaction  is  the  isolation  of  the  effect 
of  individual  species  and  obtaining  a  fundamental  knowledge  of 
specific  enhancement  processes.  Plasma  discharges  can  produce 
an  incredibly  wide  range  of  stable  and  metastable  species  that 
have  the  potential  to  enhance  combustion,  especially  when  the  dis¬ 
charge  is  in  a  premixture  of  fuel  and  oxidizer.  A  practical  system 
will  require  an  understanding  of  each  and  every  one  of  the  species 
and  their  interaction  in  a  combustion  process.  The  inclusion  of 
each  and  every  species  and  their  associated  reactions  would  render 
kinetic  mechanisms  nearly  impossible  to  use.  More  importantly, 
there  is  a  lack  of  understanding  of  what  species  and  reaction  path¬ 
ways  are  the  most  important.  Realistically,  there  will  only  be  a  few 
species  that  are  critical  to  the  enhancement  mechanism.  Therefore, 
careful  selection  of  the  plasma-produced  species  is  warranted  in 
order  to  make  the  largest  impact  on  achieving  an  understanding 
of  the  plasma-flame  interaction  and  developing  kinetic  models 
for  accurate  predictive  capabilities. 

Experiments  need  to  be  performed  to  isolate  the  individual  ef¬ 
fects  that  may  be  present  in  plasma-combustion  systems  and 
determine  to  what  extent  they  participate  in  the  enhancement 
processes.  Isolating  the  individual  enhancement  mechanisms  can 
be  difficult  because  plasma-combustion  systems  can  be  highly 
coupled  with  complex  hydrodynamics,  species  diffusion,  and  ther¬ 
mal  processes.  The  coupled  processes  do  not  allow  the  underlying 
physics  and  chemistry  involved  to  be  identified  clearly.  Further¬ 
more,  there  exists  the  additional  complication  of  plasma-produced 
species  lifetimes  when  trying  to  observe  individual  enhancement 
pathways. 

Shown  in  Fig.  1  is  an  illustration  of  some  of  the  common  plas¬ 
ma-produced  species  and  their  lifetimes  as  a  function  of  pressure. 
Depending  upon  the  application  and  the  species  of  interest,  decou¬ 
pling  the  plasma-flame  interaction  can  be  challenging.  For  exam¬ 
ple,  when  the  pressure  is  sufficiently  low  (sub-atmospheric), 
plasma-produced  species  have  long  lifetimes  and  the  potential  to 
be  transported  significant  distances.  On  the  other  hand,  at  extre¬ 
mely  high  pressures,  plasma-produced  species  have  a  very  short 
lifetime  because  of  high  rates  of  recombination  and  collisional 
quenching.  Since  the  concentration  of  species  is  proportional  to 
the  pressure  and  hence  number  density,  the  lower  the  pressure, 
the  longer  the  mean  free  path  and  the  smaller  the  collision  fre¬ 
quency.  To  minimize  the  loss  of  radicals  and  excited  species,  the 
pressures  need  to  be  low  or  the  plasma  has  to  be  located  almost  di¬ 
rectly  in  contact  with  the  reaction  zone  to  enhance  combustion. 
Unfortunately,  decoupling  of  the  plasma-combustion  interaction 
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Fig.  1.  Lifetimes  of  plasma-produced  species  as  a  function  of  pressure  in  relation  to 
common  systems  applications. 


becomes  increasingly  difficult  as  the  plasma  is  moved  closer  to 
the  reaction  zone.  Some  of  the  major  problems  that  arise  are  elec¬ 
tric  field  and  flow  disturbance  effects.  These  complications  make  it 
difficult  to  elucidate  the  fundamental  interaction,  especially  when 
trying  to  understand  the  effects  of  specific  plasma-produced 
species. 

Furthermore,  plasma  can  couple  energy  into  combustion  sys¬ 
tems  via  many  different  pathways  producing  a  diverse  range  and 
significant  number  of  species  including  radicals,  excited  species, 
ions,  and  electrons.  As  a  starting  point,  plasma-produced  species 
such  as  NOx  and  03  are  very  important,  have  been  studied  in  detail, 
and  are  integrated  with  combustion  kinetics  [5,6].  Reactive  radicals 
such  as  O,  OH,  and  H  have  been  studied  to  a  lesser  extent  in  plas¬ 
ma-assisted  combustion  systems  [19,22],  but  are  fully  integrated 
with  traditional  combustion  chemistry.  Therefore,  one  of  the  key 
links  between  combustion  and  plasma  chemistry  that  remain  are 
with  ions,  electrons,  and  excited  species.  Among  ions,  electrons 
and  excited  species,  there  has  been  some  work  on  their  effects  in 
combustion  systems,  but  the  governing  kinetics  are  weakly  under¬ 
stood  [23-25].  The  hope  is  that  the  inclusion  of  ions,  electrons,  and 
excited  species  will  open  up  the  possibility  of  new  and  faster  reac¬ 
tion  pathways  that  can  enhance  critical  combustion  phenomena 
significantly  with  a  high  degree  of  chemical  selectivity.  One  of 
the  species  that  has  attracted  much  attention  is  02(a1Ag),  also 
known  as  singlet  delta  oxygen.  It  has  a  low  excitation  energy  of 
0.98  eV  and  therefore  is  produced  in  almost  all  oxygen  containing 
plasmas.  Singlet  delta  oxygen  also  has  a  higher  oxidation  potential 
than  02  and  has  significantly  long  radiative  lifetime  (>  4000  s)  be¬ 
cause  of  a  spin  forbidden  transition  to  the  ground  state  [26]. 

There  have  been  numerous  computational  studies  aimed  to 
quantify  the  enhancement  of  ignition  and  flame  stabilization  by 
02(a1Ag).  Starik  and  Titova  showed  that  when  02  was  excited  to 
its  first  electronic  state  of  02(a1Ag)  by  laser  radiation  in  a  super¬ 
sonic  flow  of  H2-air,  the  induction  time  and  temperatures  neces¬ 
sary  for  ignition  behind  a  shockwave  were  reduced  significantly 
[27].  They  attributed  the  enhancement  to  come  from  new  path¬ 
ways  with  02(a1Ag)  to  generate  active  species,  such  as  O,  H,  and 
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OH.  Work  was  also  performed  computationally  using  an  electrical 
discharge  to  achieve  similar  results  as  the  laser  radiation  [28].  De¬ 
tailed  investigations  of  the  ignition  kinetics  with  the  presence  of 
02(a1Ag)  by  non-equilibrium  excitation  in  a  H2-02  system  was  re¬ 
viewed  in  a  paper  by  Popov  [29].  The  review  pointed  out  that  if  the 
quenching  of  02(a1  Ag)  by  H2  is  not  considered,  there  will  be  a  gross 
overestimate  of  the  amount  of  enhancement  because  the  colli- 
sional  quenching  rate  increases  significantly  with  temperature.  Po¬ 
pov  also  emphasized  was  that  there  is  a  lack  of  experimental 
studies  of  the  effect  of  electronically  excited  species  on  combustion 
phenomena. 

Starik,  Kuleshov,  and  Titova  studied  ignition  enhancement  in  an 
H2-02  mixture.  Here,  they  created  02(a1Ag)  along  with  0(1D)  with 
a  flash  photolysis  approach:  03  was  the  sensitizer  species  and  the 
radiation  was  from  a  UV  laser  beam.  The  collective  effect  of  these 
two  species  led  to  decreased  ignition  delay  times  by  several  orders 
of  magnitude  [30].  The  absorbed  energy  by  UV  laser  radiation  was 
very  small  compared  to  the  energy  required  for  the  same  amount 
of  enhancement  by  thermal  means,  indicating  the  distinct  advan¬ 
tage  of  having  strong  non-equilibrium  excitation  for  combustion 
enhancement.  Nevertheless,  the  results  were  not  verified  with 
experiments. 

Kozlov,  Starik,  and  Titova  conducted  numerical  simulations  to 
show  the  enhancement  of  H2-02  flame  speed  with  02(a1  Ag)  addi¬ 
tion  [31].  The  results  showed  that  a  concentration  of  10%  02(a1Ag) 
gave  more  than  a  50%  increase  in  the  laminar  flame  velocity.  Lean 
mixtures  were  enhanced  more  than  stoichiometric  and  rich  mix¬ 
tures  because  of  the  chain  initiation  and  branching  reactions  in¬ 
volved.  Bourig  et  al.  extended  the  numerical  modeling  by 
investigating  ignition  and  flame  propagation,  as  well  as  flame  sta¬ 
bilization  by  02(a1Ag)  [32]. 

Experimental  validation  of  the  enhancement  of  02(a1Ag)  on 
ignition  was  introduced  by  measurement  of  the  emission  from 
OH*  at  306.4  nm  [33].  Here  they  monitored  OH*  behind  a  shock- 
wave  in  an  H2-02  mixture  that  was  subjected  to  an  electrical 
discharge.  It  was  found  that  there  was  accurate  reproduction 
of  the  experimental  results  with  the  kinetic  mechanism  and  that 
OH*  might  be  a  good  indicator  for  enhancement  by  02(a1Ag)  be¬ 
cause  of  the  elevated  radical  pool  concentrations.  Skrebkov  and 
Karkach  performed  subsequent  simulations  using  the  ignition 
and  emission  spectroscopy  results  from  Ref.  [33]  and  showed 
that  there  was  reasonable  agreement;  however,  they  emphasized 
that  the  main  stumbling  block  was  the  “availability  of  experi¬ 
mentally  evaluated  amounts  of  02(a1Ag),”  as  well  as  other  plas¬ 
ma  excited  species  [34]. 

There  has  only  been  one  experimental  investigation  of  the  ef¬ 
fects  of  02(a1Ag)  on  combustion  phenomena.  Smirnov  et  al.  per¬ 
formed  experiments  aimed  at  isolating  the  effect  of  02(a1Ag)  on 
the  ignition  of  H2-02  mixtures  at  low  pressure  between  1.33  kPa 
and  2.67  kPa  [35].  The  results  showed  decreased  induction  times. 
Here,  they  noted  that  to  achieve  the  same  delay  by  heating  the 
mixture,  the  energy  input  would  be  4.4  times  greater  than  that 
needed  to  produce  02(a1Ag).  The  experiments  were  designed  to 
minimize  all  other  plasma-produced  species,  especially  O  and  03, 
with  the  injection  of  Hg  and  its  oxide  (HgO)  into  the  flow.  The  re¬ 
sults  showed  that  there  was  a  100-fold  decrease  in  the  emission 
from  0,  but  no  absolute  measurement  of  its  concentration  was  per¬ 
formed.  Furthermore,  the  03  was  not  measured  and  the  02(a1Ag) 
was  only  measured  through  emission.  Therefore,  the  measure¬ 
ments  did  not  provide  quantitative  concentrations  of  all  the  major 
species  of  enhancement  in  the  system.  Even  though  the  experi¬ 
ments  were  the  first  of  their  kind  in  an  attempt  to  isolate  the  effect 
of  02(a1Ag)  and  the  results  agreed  reasonably  well  with  their  pre¬ 
vious  calculations  [36,37]  it  was  not  certain  whether  02(a1  Ag)  was 
the  only  species  causing  the  enhancement.  In  addition,  the  exper¬ 
iments  were  conducted  only  for  an  H2-02  mixture. 


Therefore,  there  is  scarce  experimental  data  detailing  and  quan¬ 
tifying  the  effect  of  electronically  excited  species,  specifically 
02(a1Ag),  on  combustion  phenomena.  A  majority  of  the  numerical 
and  experimental  work  has  been  focused  on  the  H2-02  reaction 
system,  with  no  experimental  studies  of  the  isolated  effects  of 
02(a1Ag)  on  hydrocarbon  based  fuels.  There  have  been  no  experi¬ 
mental  studies  at  higher  pressures,  greater  than  2.67  kPa,  where 
a  structured  flame  can  exist.  Furthermore,  there  have  been  no 
experimental  studies  with  quantitative  concentration  measure¬ 
ments  of  02(a1Ag),  O,  and  03. 

The  goal  of  the  present  work  was  to  isolate  and  measure  quan¬ 
titatively  the  effects  of  02(a1Ag)  on  the  enhancement  of  flame 
propagation  using  a  hydrocarbon  based  fuel,  as  well  as  to  under¬ 
stand  the  kinetic  mechanisms  involved.  The  enhancement  of  flame 
propagation  speeds  were  investigated  through  the  development  of 
an  integrated  plasma-combustion  experimental  platform  where 
02(a1Ag)  was  produced,  transported,  and  absolute  measurements 
taken  through  integrated-cavity-output  absorption  spectroscopy. 
A  low  power  microwave  coupled  plasma  was  used  to  activate  Ar/ 
02  as  the  oxidizer  co-flow  of  low  pressure  C2H4  lifted  flames.  The 
individual  effects  of  02(a1Ag)  were  isolated  by  NO  addition  to  cat- 
alytically  remove  03  and  O,  the  two  major  quenching  species  of 
02(a1Ag).  The  experiments  provide  the  first  experimental  evidence 
of  the  isolated  effect  of  02(a1Ag)  on  the  propagation  of  a  hydrocar¬ 
bon  fuel  based  flame.  The  results  will  provide  important  steps  to¬ 
wards  developing  a  comprehensive  predictive  model  for  plasma- 
combustion  systems  with  detailed  and  well-defined  experimental 
results. 

2.  Experimental  system 

The  experimental  platform  was  designed  to  have  two  separate 
apparatuses  consisting  of  a  lifted  flame  burner  and  a  plasma-pro¬ 
duced  species  diagnostics  system.  The  two  apparatuses  shared 
the  same  flow  controllers,  plasma  device,  and  inert  wall  surfaces. 
A  schematic  of  the  experimental  platform  is  shown  in  Fig.  2. 

2.1.  Plasma-assisted  lifted  flame  system 

A  laminar  lifted  flame  burner  was  adopted  for  the  combustion 
platform  and  was  placed  in  a  variable  pressure  chamber  that  could 
be  used  from  atmospheric  pressure  down  to  2.67  kPa  with  the  in¬ 
stalled  vacuum  and  flow  system.  The  lifted  flame  burner  consisted 
of  a  central  fuel  jet  with  a  diameter  of  1.04  mm  and  was  located  in 
a  90  mm  inner  diameter  fused  silica  (quartz)  tube  to  contain  the 
co-flow  of  02  and  Ar  in  the  chamber.  The  fuel  nozzle  was  aerody- 
namically  shaped  to  produce  a  uniform  velocity  profile  at  the  exit. 
To  ensure  that  the  co-flow  was  uniform,  two  stainless  steel  meshes 
coated  with  quartz  for  chemical  inertness  were  separated  by  3  cm 
and  were  located  between  the  oxidizer  inlet  of  the  burner  and  the 
fuel  jet  exit.  The  gases  used  in  the  experiments  were  C2H4  for  the 
fuel  and  ultra-high  purity  02  (99.99%)  and  Ar  (99.95%)  mixed  for 
the  oxidizer.  The  flow  rate  of  the  fuel  was  controlled  with  a  cali¬ 
brated  mass  flow  meter  while  the  02  and  Ar  were  controlled  with 
calibrated  sonic  nozzles  to  give  flow  rate  uncertainties  of  less  than 
1%. 

To  excite  the  co-flow  in  the  low  pressure  experiments,  an  elec¬ 
trodeless  microwave  discharge  (McCarroll  cavity  driven  by  an 
Opthos  MPG-4M  microwave  power  supply)  with  up  to  100W  of 
power  was  used  external  to  the  chamber,  upstream  of  the  lifted 
flame  burner  to  activate  the  02  in  the  mainly  Ar  flow  (15%  02  in 
85%  Ar  for  3.61  kPa  and  11.9%  02  in  88.1%  Ar  for  6.73  kPa).  The 
plasma  was  initiated  in  the  microwave  cavity  by  seeding  the  up¬ 
stream  flow  with  ionized  gas  created  by  a  high-voltage,  high-fre¬ 
quency  Tesla  coil.  The  mixture  of  ultra-high  purity  02/Ar  was 
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Fig.  2.  Experimental  set-up  schematic  of  variable  pressure  lifted  flame  burner  integrated  with  microwave  plasma  discharge  device  and  plasma  flow  diagnostics  system.  The 
abbreviation  “PD”  denotes  photodetector. 


therefore  activated  by  the  self-sustained  microwave  discharge 
which  was  maintained  when  the  Tesla  coil  was  switched  off.  The 
plasma  system  was  chosen  because  of  its  flexibility  of  being  used 
external  to  a  quartz  tube  flow  system,  as  well  as  its  ease  of  tuning 
and  stability  for  the  range  of  pressure  used  in  the  experiments.  Fur¬ 
thermore,  the  lower  power  output  of  the  plasma  system  produced 
a  glow  discharge  with  lower  reduced  electric  field  and  an  electron 
energy  distribution  function  that  was  peaked  at  lower  electron 
temperature  where  significant  concentrations  of  02(a1Ag)  could 
be  excited  from  ground  state  02.  Nevertheless,  the  microwave  dis¬ 
charge  produced  excited  Ar,  as  well  as  multiple  oxygen  containing 
species  including  O,  03,  02(v),  0(1D),  0(1S),  02(a1Ag),  O^b1^),  etc. 

The  high  velocity  fuel  jet  (approximately  20-40  m/s)  and  low 
velocity  co-flow  (approximately  0.15-0.2  m/s)  created  a  mixing 
layer  with  a  stoichiometric  contour  where  the  premixed  flame 
head  of  a  lifted  flame  was  located  (shown  in  the  top  right  inset 
in  Fig.  2).  The  lifted  flame,  which  is  also  called  a  tribrachial  (triple) 
flame,  had  a  premixed  flame  head  and  diffusion  flame  tail  where 
the  flame  was  always  anchored  on  the  stoichiometric  contour. 


The  lifted  flame  could  be  located  at  different  distances  from  the 
fuel  jet  nozzle  depending  upon  the  local  flow  velocity.  For  a  fixed 
flow  field,  the  flame  is  located  in  a  stationary  position  where  the 
lifted  flame  speed  at  the  premixed  flame  head  is  balanced  with 
the  local  flow  velocity.  In  Fig.  3,  photographs  are  shown  of  station¬ 
ary  C2H4  lifted  flames  as  a  function  of  fuel  jet  velocity  for  a  pressure 


Fig.  3.  Photographs  of  C2H4  lifted  flames  at  stationary  positions  for  different  fuel  jet 
velocities  at  6.73  kPa. 
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of  6.73  kPa.  If  the  flame  speed  increases,  the  liftoff  height  decreases 
to  re-establish  a  local  dynamic  balance  between  the  flame  speed 
and  flow  velocity. 

Due  to  the  slow  laminar  boundary  layer  development  and  the 
velocity  and  concentration  gradients  created,  the  lifted  flame 
height  is  very  sensitive  to  the  changes  in  flame  speed  and  therefore 
provides  excellent  flame  geometry  for  the  direct  observation  of 
flame  speed  enhancement.  Since  the  fuel  and  oxidizer  are  not 
mixed  far  upstream  of  the  flame,  there  is  very  short  residence  time 
for  the  fuel  and  oxidizer  to  react  in  the  cold  flow.  The  short  resi¬ 
dence  time  helps  to  further  decouple  the  enhancement  effects  to 
be  directly  from  reactions  in  the  flame  zone  and  not  far  upstream 
in  the  cold  un-reacted  flow. 

All  surfaces  that  the  plasma  afterglow  gases  came  in  contact 
with  were  treated  to  be  chemically  inert.  The  tubing  used  was 
fused  silica  (quartz),  the  fittings  were  316  stainless  steel  (non¬ 
magnetic)  and  the  lifted  flame  burner  was  coated  with  silica  (Res- 
tek  Silcosteel).  The  inert  surfaces  mitigated  the  quenching  from  ac¬ 
tive  wall  surfaces  and  promoted  the  transport  of  02(a1Ag)  to  the 
flame.  Multiple  temperatures  and  pressures  were  monitored  and 
recorded  in  the  system  with  thermocouples  and  pressure  transduc¬ 
ers,  respectively.  Temperatures  were  measured  at  points  Ti,  T2,  and 
T3  shown  in  Fig.  2  corresponding  to  the  burner  inlet  tube  surface 
temperature,  burner  surface  temperature,  and  co-flow  gas  temper¬ 
ature,  respectively.  Additionally,  the  pressure  upstream  of  the 
microwave  plasma  and  in  the  chamber  respectively  at  points  Pi 
and  P2  in  Fig.  2  were  continuously  monitored  and  recorded. 

2.2.  Plasma  afterglow  diagnostics  systems 

The  plasma  activated  oxidizer  flow  was  run  through  a  series  of 
diagnostics  to  quantify  the  concentrations  of  species  produced  by 
the  plasma.  The  flow  system  used  the  same  chemically  inert  flow 
surfaces  and  residence  times  for  a  direct  comparison  to  the  lifted 
flame  system.  The  diagnostics  were  not  run  in  situ  with  the  flame 
system,  but  instead  as  two  separate  systems  with  common  flow 
control  and  plasma  discharge.  The  systems  were  used  separately 
because  having  the  ICOS  cavity  between  the  plasma  and  flame 
would  double  the  transport  residence  time  and  decrease  the  con¬ 
centrations  of  the  species  of  interest.  For  the  3.61  kPa  and 
6.73  kPa  experiments,  the  residence  time  from  the  plasma  to  the 
end  of  the  ICOS  cavity  was  approximately  1320  ms  and  1680  ms, 
respectively  and  between  the  plasma  and  flame  in  the  burner  were 
approximately  910  ms  and  1470  ms,  respectively. 

2.2.1.  Quantitative  02(a1Ag)  Measurement  Using  Off-Axis  ICOS 

The  02(a1  Ag)  produced  by  the  microwave  plasma  discharge  was 
measured  by  using  highly  sensitive  integrated-cavity-output  spec¬ 
troscopy  by  absorption  at  the  (1,0)  band  of  the  b1^/  -  a1Ag  Nox- 
on  system  [38].  The  ICOS  system  measured  the  average  number 
density  of  02(a1Ag)  across  an  82.5  cm  long  (approximately 
1500  ms  residence  time)  absorption  cell  downstream  of  the  plas¬ 
ma.  The  effective  path  length  was  greater  than  78  km  due  to  multi¬ 
ple  passes  and  provided  accurate  measurements  down  to 
1014  molecules/cm3.  In  Fig.  4  the  transmittance  spectrums  from 
the  ICOS  cavity  for  the  experimental  conditions  with  and  without 
NO  addition  and  plasma  activation  are  shown.  The  absorption  fea¬ 
ture  demarcated  by  the  red  box  shows  the  location  of  the  Q(12) 
transition  of  02(a1  Ag),  which  was  the  primary  transition  of  interest 
[38].  The  Q(12)  transition  was  chosen  for  absorption  measure¬ 
ments  because  there  was  no  interference  between  it  and  any  of 
the  absorption  features  of  other  species  present  in  the  flow.  In 
Fig.  5,  the  measured  absorption  profile  of  02(a1Ag)  using  the 
Q(12)  transition  for  plasma-activated  Ar/02  at  3.61  kPa  is  shown. 
The  fitted  curve  in  Fig.  5  was  a  function  of  the  absorption  path- 
length  and  cross-section,  as  well  as  broadening  from  pressure 


Fig.  4.  ICOS  transmittance  spectrum  with  and  without  NO  addition  and  plasma 
activation  (P  and  NP  are  plasma  and  no  plasma,  respectively). 
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Fig.  5.  Cavity-enhanced  absorption  (GA)  measurement  of  02(a1Ag)  obtained  by 
tuning  the  laser  over  the  Q(12)  transition. 

and  temperature.  A  more  detailed  description  of  the  measurement 
and  curve  fitting  process  to  obtain  absolute  number  densities  of 
02(a1Ag)  can  be  found  in  previous  work  [38]. 

2.2.2.  Quantitative  absorption  measurements  of  03 

The  production  of  02(a1Ag)  also  resulted  in  significant  produc¬ 
tion  of  03,  and  therefore  it  had  to  be  measured  quantitatively.  To 
accomplish  this,  a  one-pass,  line-of-sight  absorption  cell  was 
adopted.  The  absorption  cell  was  comprised  of  a  stainless  steel 
compression  cross  fitting  with  the  side  arms  made  of  quartz  tubes 
capped  with  UV  quality  windows.  At  one  window,  a  mercury  light 
with  stable  output  provided  ultraviolet  light  at  the  wavelengths 
covering  the  Hartley  band  of  03.  A  10  nm  notch  filter  isolating  only 
the  253.7  nm  mercury  line  was  placed  at  the  exit  of  the  absorption 
cell,  followed  by  a  photodiode  detector.  The  wavelength  of 
253.7  nm  was  used  because  03  has  a  peak  absorption  cross-section 
there  of  1.137  x  10“17  cm2  (at  300  K)  [39].  No  other  species  pres¬ 
ent  in  the  flow  absorb  at  this  wavelength.  Therefore,  the  change 
in  the  transmittance  of  the  cell  with  the  plasma  on  and  off  could 
be  used  to  determine  the  03  concentration  through  the  Beer-Lam- 
bert  law 


where  Nozone  is  the  absolute  number  density  of  the  absorbing  spe¬ 
cies,  03,  /  the  intensity  of  light  with  the  presence  of  03,  J0  the  inten¬ 
sity  of  light  without  the  presence  of  03,  (jozone  the  absorption  cross- 
section  of  03  at  the  excitation  wavelength  (253.7  nm),  and  L  the 
path  length  in  the  absorption  cell  (12.48  cm).  The  concentration  was 
then  calculated  in  parts  per  million  (ppm)  of  03  with  an  uncertainty 
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of  approximately  ±2%  and  a  minimum  detectable  threshold  of 
approximately  15  ppm.  The  uncertainty  and  minimum  detectable 
threshold  came  from  the  fluctuations  in  the  intensity  of  the 
mercury  light  as  a  function  of  time  during  the  experiments. 

2.2.3.  FTIR  spectroscopy  for  NO  and  N02  measurement 

The  isolation  of  02(a1Ag)  in  the  plasma  afterglow  was  achieved 
with  the  addition  of  a  small  concentration  of  NO  to  the  flow.  The 
details  of  this  approach  are  described  in  Section  2.3.  Therefore, 
quantitative  measurements  of  NO  and  N02  were  required  to  con¬ 
firm  the  presence  of  these  species  in  the  system.  To  accomplish 
these  measurements  the  flow  was  split  downstream  of  the  03 
absorption  cell  to  allow  sampling  with  Fourier  Transform  Infrared 
(FTIR)  spectroscopy  (Fig.  2).  The  pressure  and  temperature  were 
fixed  at  2  kPa  and  373  K,  respectively,  in  the  FTIR  cell  for  all  exper¬ 
imental  measurements  to  maintain  sufficient  flow  rates  and  sam¬ 
pling  times.  The  absorption  features  of  NO  and  N02  were  chosen 
where  there  would  be  no  chance  of  interference  from  changes  in 
the  background  or  other  species. 

2.3.  Isolation  of  02(alAg)  and  03from  other  plasma-produced  species 

To  avoid  the  complications  of  plasma  disturbances  in  the  com¬ 
bustion  system,  the  flow  was  activated  far  upstream  of  the  lifted 
flame.  With  the  flow  rates  used  in  the  experiments,  the  average 
residence  times  between  the  plasma  and  the  flame  was  approxi¬ 
mately  1  s.  The  significantly  long  residence  time  was  chosen  in  or¬ 
der  to  quench  the  plasma-produced  species  that  were  not  of 
interest  in  the  experiments.  Furthermore,  since  the  focus  of  the 
current  experiments  were  to  isolate  02(a1Ag),  Ar  was  used  as  the 
inert  instead  of  N2.  Using  Ar  served  two  purposes.  First,  Ar  has  only 
electronic  excitation  requiring  11.6  eV  for  the  first  electronic  level. 
Therefore,  when  using  the  low  power  microwave  system,  most  of 
the  energy  will  then  be  deposited  into  the  02  to  produce  O,  03, 
02(a1Ag),  02(b1Sg),  and  metastable  O  and  Ar.  Second,  without  N2, 
there  would  be  no  nitrogen  containing  species  produced,  specifi¬ 
cally  NO  or  N02,  simplifying  the  chemistry  in  the  plasma  and  in 
the  post-plasma  flow. 

With  the  1  s  residence  time,  the  only  two  species  that  have  long 
enough  lifetime  to  allow  for  measurement  and  introduction  to  a 
combustion  system  are  03  and  02(a1Ag).  A  list  of  the  primary 
quenching  reactions  of  the  oxygen  containing  species  is  given  in 
Table  1  in  order  of  reaction  rate.  Quenching  reactions  with  02  dom¬ 
inate  over  the  inert  species  present  and  are  therefore  the  reactions 
listed  for  a  comparison.  Ozone  is  the  most  stable  with  02(a1  Ag)  also 
having  long  lifetime.  All  other  species  are  quenched  many  times 
faster. 

Beyond  the  gas  phase  kinetics,  the  wall  quenching  effects  have 
to  be  considered.  Many  species  can  be  transported  significant  dis¬ 
tances  with  inert  wall  surfaces.  Therefore,  by  coating  all  wall  sur¬ 
faces  with  quartz  or  using  quartz  tubes  is  essential  to  ensure 
minimal  quenching.  Some  examples  of  the  relative  reaction  prob¬ 
abilities  of  plasma-produced  species  with  a  wall  surface  of  Pyrex 
(similar  to  quartz)  are  shown  in  Table  2.  Therefore,  having  large 


Table  1 

Reaction  rates  of  plasma-produced  oxygen  species  at  298  K.  The  term  “HP”  refers  to 
the  high  pressure  limit. 


Reaction 

Reaction  constant  (cm3/molecule/s) 

0  +  02  +  M->03  +  M 

6.0  x  10“34  =  (HP  limit  3.61  x  10”10)  [40] 

0(’D)  +  02^0  +  02 

4.0  x  10“n  [40] 

02(v)  +  02  — ►  02  +  02 

1.73  x  10“13  [41] 

02(b1Eg)  +  02  ->  02  +  02 

4.1  x  10“17  [40] 

02(a1Ag)  +  02  ->  02  +  02 

1.6  x  10“18  [40] 

O^a1  Ag)  +  Ar  -►  02  +  Ar 

1.0  x  IQ”20  [42] 

Table  2 

Wall  quenching  probabilities  of  plasma-produced  species  for  a  Pyrex  surface. 


Wall  reaction 

Reaction  probability 

02(a1Ag)  +  wall  -»•  02 

2  x  10“5 [43] 

02(b1Zg)  +  wall  ->  02 

2  x  10“2 [43] 

02(v)  +  wall  — ►  02 

0.2  [41] 

0(1D)  +  wall^O(3P) 

1.0  [43] 

0(1S)  +  wall  ->  0(3P) 

1.0  [43] 

O  +  wall  ->  Vi02 

2  x  10“2 [43] 

M+  +  wall  -»•  M 

1.0  [41] 

surface  to  volume  ratios  in  the  flow  will  help  to  suppress  the  con¬ 
centrations  of  all  species  except  03  and  02(a1Ag). 

There  remained  the  need  to  isolate  03  and  02(a1Ag)  to  observe 
their  individual  effects.  Ozone  can  be  isolated  by  simply  using 
higher  pressures  because  the  02(a1Ag)  will  be  collisionally 
quenched  to  02(X3Zg).  This  was  demonstrated  in  part  I  of  this 
two  part  work  [44].  To  isolate  02(a1  Ag),  a  different  approach  needs 
to  be  taken  to  remove  03.  To  mitigate  the  effect  of  03  and  com¬ 
pletely  isolate  the  effect  of  02(a1Ag),  NO  was  added  to  the  flow 
downstream  of  the  plasma  in  the  flow  diagnostics  system,  as  well 
as  the  combustion  system.  The  addition  of  NO  in  prescribed  con¬ 
centrations  served  two  purposes.  First,  it  catalytically  removed 
03  from  the  system  to  isolate  02(a1Ag)  and  second,  NO  addition 
in  the  concentrations  in  the  experiments  would  not  interfere  with 
the  effects  of  02(a1Ag). 

The  isolation  of  02(a1Ag)  with  NO  addition  relies  upon  the  reac¬ 
tion  of  NO  with  03  via 

NO  +  03  N02  +  02  (2) 

Reaction  (2)  is  over  three  orders  of  magnitude  faster  than  NO  with 
02(a1Ag),  as  shown  in  Table  3.  Furthermore,  the  subsequent  reac¬ 
tion  of  N02  with  02(a1Ag)  is  also  slow  in  comparison  to  reaction 
(2).  The  major  consumption  pathway  of  N02  would  be  from  the 
reaction  with  0  atoms  via 

N02  +  0  NO  +  02  (3) 

Since  0  and  03  are  present  at  the  same  location  in  the  flow  tube,  i.e., 
0  atoms  are  converted  to  03,  NO  acts  as  a  catalyst  and  very  little  is 
needed  to  fully  eliminate  the  03.  Therefore,  the  individual  effects  of 
02(a1Ag)  can  be  selected  by  the  use  of  NO  addition.  Finally,  the  pres¬ 
ence  of  NO  can  also  be  used  to  determine  the  presence  of  0  atoms 
via  the  well  known  0  atom  titration  reaction 

NO  +  0  ->  NO;  N02  +  hv  (4) 

The  absence  of  the  greenish-yellow  glow  in  the  plasma  after¬ 
glow  is  used  to  assure  that  all  the  O  atoms  are  quenched  and  not 
present  in  the  flow  system. 

3.  Results  and  discussion 

3A.  NO  Addition  for  02(a^Ag)  Isolation 

To  investigate  the  flame  speed  enhancement  by  02(a1Ag),  re¬ 
duced  pressures  were  used  to  suppress  the  quenching  and  recom- 


Table  3 

Reaction  rates  of  02(a1Ag)  and  03  with  NO  and  N02  at  298  K. 


Reaction 

Reaction  constant  (cm3/molecule/s) 

O^Ag)  +  NO  ->  02  +  NO 

4.48  x  10“17  [45] 

02(a1Ag)  +  NO  ->  O  +  N02 

4.88  x  10“18  [42] 

Oz^Ag)  +  N02  ->  02  +  N02 

5.00  x  10-18  [42] 

03  +  NO  ->  02  +  N02 

1.80  x  10“14  [40] 

03  +  N02  — ►  02  +  02  +  NO 

1.00  x  10“18  [46] 
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Fig.  6.  Experimentally  measured  concentrations  of  03  and  02(a1Ag)  as  a  function  of 
pressure  before  03  reduction  by  NO  addition. 

bination  rates.  The  ICOS  and  03  absorption  diagnostics  were  used 
to  initially  measure  02(a1Ag)  and  03  as  a  function  of  pressure 
and  02  concentrations.  The  results  in  Fig.  6  show  that  there  remain 
significant  concentrations  of  03,  despite  the  decrease  of  the  pres¬ 
sure.  At  lower  pressure  and  02  loading  in  Ar,  the  02(a1Ag)  was  in 
significantly  greater  concentration  than  03.  Unfortunately,  stable 
flames  under  lower  pressure  and  02  loadings  were  difficult  to 
achieve.  Therefore,  the  removal  of  03  was  warranted  for  the  iso¬ 
lated  study  of  02(a1Ag). 

The  results  in  Fig.  7  show  that  the  addition  of  NO  suppressed 
the  concentration  of  03  to  a  magnitude  below  the  threshold  that 
could  be  measured  (ppm  levels)  while  the  02(a1Ag)  concentration 
at  a  given  residence  time  increased  by  almost  an  order  of  magni¬ 
tude  over  a  wide  pressure  range.  Since  the  NO  addition  worked  cat- 
alytically  to  reduce  the  02(a1Ag)  quenching  species  of  03  and  O 
through  reactions  (2)  and  (3),  the  process  was  heavily  reliant  upon 
the  concentration  of  O  where  the  NO  was  added  to  the  system. 
With  the  flow  rates  used  in  the  system,  the  residence  time  between 
the  plasma  and  the  location  of  NO  addition  was  approximately 
5  ms.  For  the  pressures  used  in  the  experiments,  there  was  a  signif¬ 
icant  concentration  of  O  in  relation  to  03.  Nevertheless,  the  concen¬ 
trations  of  NO  and  N02  needed  to  be  verified  to  ensure  that  the  N02 
concentration  remained  negligible. 

To  confirm  that  there  was  no  conversion  of  NO  to  N02,  measure¬ 
ments  were  taken  using  the  FTIR  with  the  plasma  off  and  on.  The 
measurements  were  taken  initially  of  the  plasma  off  where  the 
FTIR  measured  the  exact  concentration  of  NO  that  was  being  added 
to  the  system.  When  the  plasma  was  turned  on,  the  NO  concentra¬ 
tion  did  not  change,  as  shown  in  Fig.  8.  To  confirm  further  that 


Fig.  7.  Experimentally  measured  concentrations  of  02(a1Ag)  and  03  with  and 
without  500  ppm  NO  for  15%  02  in  85%  Ar  as  a  function  of  pressure. 


Frequency  [cm  J] 

Fig.  8.  NO  absorption  measurements  from  FTIR  showing  that  the  NO  concentra¬ 
tions  did  not  change  with  the  plasma  on  or  off  and  therefore  no  conversion  to  N02 
(P  and  NP  are  plasma  and  no  plasma,  respectively). 


2920.5  2921.2  2921.9  2922.6  2923.3  2924 

Frequency  [cm1] 

Fig.  9.  N02  absorption  measurements  from  FTIR  showing  that  the  NO  is  not 
converted  to  N02  when  the  plasma  is  on  or  off  (P  and  NP  are  plasma  and  no  plasma, 
respectively). 


there  was  no  N02  production,  Fig.  9  shows  that  the  concentration 
of  N02  remained  below  the  detectability  threshold  of  the  FTIR. 
Two  reference  absorption  spectrums  of  500  ppm  and  4000  ppm 
of  N02  that  were  taken  at  the  same  pressures  and  temperatures 
are  shown  for  comparison  and  that  the  three  peaks  present  can 
be  resolved.  The  results  show  clearly  that  the  highest  concentra¬ 
tion  of  N02  was  on  the  order  of  10’s  of  ppm.  Since  flame  speed 
enhancement  by  N02  is  more  than  a  factor  of  three  smaller  than 
the  effect  of  03  for  the  same  concentrations,  the  small  concentra¬ 
tion  of  N02  would  have  a  negligible  effect  on  flame  speed.  There¬ 
fore,  when  NO  was  added  to  the  plasma  afterglow,  the  only 
plasma-produced  species  to  survive  would  be  02(a1Ag)  and  the 
NO  concentration  would  remain  constant. 

To  confirm  further  the  chemical  kinetic  processes  involved  in 
the  post-discharge  gases  in  the  system,  a  simple  flow  kinetic  model 
was  compiled  to  include  O,  03,  02(a1  Ag),  N2,  NO,  N02,  and  Ar.  Other 
excited  species  were  excluded  because  they  existed  in  only  small 
concentrations  at  the  point  of  NO  injection,  which  was  the  starting 
point  for  the  numerical  simulations.  Specifically,  02(b1Zg)  was  ex¬ 
cluded  because  it  would  be  more  than  an  order  of  magnitude  lower 
in  concentration  than  02(a1Ag)  and  would  also  quench  rapidly 
[35].  Nitrogen  was  included  in  the  model  because  the  NO  being  in¬ 
jected  into  the  system  downstream  of  the  plasma  was  accompa¬ 
nied  by  N2.  The  presence  of  N2  did  not  affect  any  of  the  results 
because  it  did  not  react  with  any  of  the  other  species  present  in 
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the  system.  The  most  critical  reaction  would  be  with  02(a1Ag),  but 
the  quenching  rate  is  of  the  same  order  as  the  rate  with  Ar  [43]. 
One  of  the  most  important  initial  conditions  for  the  numerical  sim¬ 
ulations  was  the  concentration  of  atomic  oxygen.  To  find  what  con¬ 
centrations  existed  in  the  experiments,  an  N02  titration  technique 
was  adopted  [47].  The  N02  was  added  to  the  plasma  afterglow  at 
the  same  location  of  NO  injection,  approximately  5  ms  down¬ 
stream  of  the  plasma.  The  N02  titration  technique  works  by  reac¬ 
tion  (3)  being  five  orders  of  magnitude  faster  than  reaction  (4). 
Therefore,  when  N02  was  injected  after  the  plasma,  the  FTIR  was 
used  to  sample  the  flow  downstream  and  monitor  the  NO  versus 
N02  concentration.  The  N02  was  continually  added  until  the  FTIR 
showed  the  presence  of  N02  and  no  changes  in  the  concentration 
of  NO.  At  that  point,  the  concentration  of  N02  being  added  was 
equal  to  the  concentration  of  O  at  the  injection  location  in  the  flow. 
Furthermore,  the  lack  of  emission  from  reaction  (4)  would  also 
indicate  that  there  was  no  longer  O  present  in  the  flow.  The  N02 
titration  technique  provided  quantification  of  the  O  concentration 
within  a  10%  uncertainty  to  be  used  as  an  initial  condition  in  the 
kinetic  model. 

Without  NO  addition  to  the  plasma  afterglow,  the  model  pre¬ 
dicted  the  concentrations  of  03  and  02(a1Ag)  well.  With  NO  addi¬ 
tion,  it  was  found  that  if  the  simple  flow  kinetic  model  was 
maintained  at  300  K,  the  catalytic  cycle  to  remove  03  was  not  com¬ 
plete  and  produced  significant  concentrations  of  N02.  The  results 
indicated  that  reaction  (3)  was  not  participating,  possibly  from 
the  lack  of  presence  of  O.  Realizing  that  the  temperature  of  the 
gas  in  the  plasma  afterglow  is  not  300  K,  but  starts  at  a  tempera¬ 
ture  around  450  K  and  decreases  to  nearly  room  temperature  after 
a  short  residence  time,  the  model  was  adjusted  to  account  for  this. 
With  a  prescribed  temperature  gradient,  the  results  with  NO  addi¬ 
tion  showed  clearly  that  the  catalytic  cycle  did  not  consume  any 
NO  and  the  N02  remained  negligible  (Fig.  10).  Furthermore,  the  re¬ 
sults  without  NO  addition  agreed  well  with  the  experiments 
(Fig.  10).  The  temperature  gradient  was  critical  to  the  catalytic  cy¬ 
cle  because  of  the  quenching  reactions  of  atomic  oxygen.  When  the 
temperature  was  fixed  at  300  K,  the  O  quenched  quickly  to  produce 
03,  therefore  not  allowing  the  N02  to  be  converted  back  to  NO. 
When  the  temperatures  were  higher  at  the  beginning  of  the  com¬ 
putation,  the  recombination  reactions  were  suppressed,  allowing 
for  the  catalytic  cycle  to  complete. 

The  validation  of  the  simple  flow  model  with  the  experimental 
results  allowed  for  it  to  be  used  to  find  the  concentrations  of 
02(a1Ag)  at  the  residence  time  of  where  the  flame  was  in  the  bur¬ 
ner  system.  Furthermore,  the  change  of  02(a1Ag)  concentration  at 


Fig.  10.  Plot  of  numerical  simulation  results  of  plasma  afterglow  species  profiles 
with  and  without  NO  addition. 


the  residence  time  of  the  flame  was  very  small,  allowing  for  mini¬ 
mal  error.  Therefore,  the  02(a1Ag)  could  be  measured  in  the  ICOS 
cavity  and  the  kinetic  model  used  to  find  the  actual  concentration 
at  the  specific  flow  residence  time  where  the  flame  was  located  in 
the  combustion  system. 

3.2.  Experimental  results  of  flame  propagation  enhancement  by 
02(a'Ag) 

To  examine  the  effects  of  02(a'Ag)  on  flame  propagation  speed, 
C2H4  was  used  as  the  fuel  to  produce  stable  lifted  flames  at  low 
pressures.  The  co-flow  conditions  of  velocity  and  02  concentration 
in  Ar  were  fixed,  along  with  the  fuel  jet  velocity  to  establish  a  flame 
at  a  stationary  lifted  location.  Photographs  were  taken  of  the  flame 
while  simultaneously  recording  the  pressure  and  temperature.  The 
microwave  plasma  was  then  turned  on  and  photographs  taken 
again  of  the  flame  with  a  lower  liftoff  height.  The  change  in  flame 
liftoff  height  was  calculated  from  the  photographs.  A  similar  proce¬ 
dure  was  used  again  with  NO  addition  just  downstream  of  the 
microwave  plasma  cavity,  replicating  the  residence  times  in  the 
plasma  afterglow  diagnostics  system.  The  results  showed  that 
there  was  a  significant  change  in  the  flame  liftoff  height  when 
the  plasma  was  turned  on  with  and  without  the  addition  of  NO. 
Therefore,  the  flame  propagation  enhancement  came  from  a  com¬ 
bination  of  03  and  02(a1Ag)  without  NO  addition  and  only  from 
02(a1Ag)  with  NO  addition.  Concentrations  between  500  ppm 
and  2000  ppm  of  NO  were  added  downstream  of  the  plasma  to  give 
different  concentrations  of  02(a1Ag)  at  the  flame  front.  The  more 
NO  that  was  added,  the  faster  the  03  and  O  were  removed  before 
they  reacted  with  and  quenched  02(a1Ag). 

The  experiments  were  performed  for  several  conditions  at  both 
3.61  kPa  and  6.73  kPa.  The  flow  field,  temperature,  and  pressure  re¬ 
mained  constant  when  the  plasma  was  cycled  off  and  on,  so  a  di¬ 
rect  comparison  between  plasma  off  and  on  can  be  attributed  to 
the  enhancement  by  03  and  02(a1Ag).  The  temperature  and  pres¬ 
sure  remained  constant  within  an  uncertainty  of  0.1  K  and 
26.7  Pa  respectively.  Experiments  were  performed  to  quantify  that 
the  uncertainty  in  temperature  and  pressure  was  not  affecting  the 
flame  liftoff  height  enough  to  mask  the  enhancement  by  03  and 
02(a1  Ag)  addition.  For  a  temperature  change  of  0.1  K,  the  flame  lift¬ 
off  height  changed  by  0.07  mm  and  for  a  pressure  change  of  26.7  Pa 
the  flame  liftoff  height  changed  by  0.29  mm.  Therefore,  for  the 
change  in  flame  liftoff  height  observed  in  the  experiments  with 
03  and  02(a1Ag)  addition,  which  was  on  the  order  of  5  mm  to 
10  mm,  the  uncertainty  was  more  than  an  order  of  magnitude 
smaller. 

Also  of  significance  to  note  was  that  the  presence  of  NO  did  af¬ 
fect  the  flame  structure  by  changing  the  stoichiometric  contour 
and  flame  speed  slightly,  but  the  comparison  was  between  the 
plasma  being  off  and  on  with  constant  NO  addition.  Furthermore, 
since  the  flow  diagnostics  showed  that  the  NO  acted  catalytically 
to  reduce  the  03  and  did  not  change  in  concentration,  the  condi¬ 
tions  of  the  plasma  being  off  and  on  could  be  compared  directly 
to  each  other.  The  near  zero  concentrations  of  N02  would  have  a 
negligible  effect  on  flame  speed  since  the  enhancement  was  calcu¬ 
lated  to  be  more  than  three  times  lower  than  that  of  03. 


Table  4 

Change  in  flame  liftoff  height  (A HL)  with  simulation  corrected  concentrations  of 
02(a1  Ag)  and  03  present  at  flame  for  a  plasma  power  of  80  W  at  3.61  kPa. 


02(a1Ag)  (ppm) 

03  (ppm) 

A HL  (mm) 

215 

513 

7.15 

3851 

0 

4.76 

5416 

0 

7.31 

5571 

0 

6.82 

5596 

0 

6.83 
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Table  5 

Change  in  flame  liftoff  height  (A HL)  with  simulation  corrected  concentrations  of 
O^a1  Ag)  and  03  present  at  flame  for  a  plasma  power  of  80  W  at  6.73  kPa. 


02(a1Ag)  (ppm) 

03  (ppm) 

A Hh  (mm) 

10 

423 

7.70 

2345 

0 

1.91 

3285 

0 

5.64 

3391 

0 

6.11 

The  same  02  loadings,  flow  rates,  and  pressures  were  then  used 
on  the  plasma  afterglow  diagnostics  system  to  find  the  concentra¬ 
tions  of  03  and  02(a1Ag)  that  were  present  at  the  lifted  flame.  The 
measured  concentrations  of  03  and  02(a1Ag)  were  corrected  for 
residence  time  using  the  kinetic  model  described  in  Section  3.1. 
The  results  are  shown  in  Tables  4  and  5  as  a  function  of  03  and 
02(a1Ag)  concentration  for  3.61  kPa  and  6.73  kPa,  respectively. 
Overall  with  NO  addition,  there  was  no  change  in  the  flow  field, 
temperature,  or  any  species  other  than  02(a1Ag)  when  the  plasma 
was  cycled  on  and  off.  Therefore,  by  turning  on  the  plasma,  it  was 
the  equivalent  of  introducing  a  pure  source  of  02(a1  Ag).  The  results 
in  Tables  4  and  5  show  that  there  was  a  clear  correlation  between 
the  concentration  of  02(a1Ag)  present  and  the  change  in  flame  lift¬ 
off  height.  Fig.  1 1  shows  the  trend  of  increasing  change  in  flame 
liftoff  height,  and  hence  more  flame  speed  enhancement,  with 
increasing  02(a1Ag)  concentration  for  both  3.61  kPa  and  6.73  kPa. 
The  experimental  uncertainties  for  the  change  in  flame  liftoff 
height  was  ±0.5  mm  and  for  02(a1Ag)  concentration  was 
±500  ppm,  and  are  shown  by  the  error  bars  in  Fig.  11. 

With  the  results  showing  clearly  the  enhancement  of  flame 
propagation  speed  by  02(a1Ag),  it  was  important  to  determine 
the  enhancement  quantitatively.  Unlike  the  lifted  flame  at  atmo¬ 
spheric  pressure  described  in  the  first  part  of  this  work  related  to 
03,  a  cold  flow  similarity  solution  does  not  correctly  describe  the 
low  pressure  experiments.  Therefore,  an  indirect  method  was  used 
to  find  the  amount  of  flame  speed  enhancement  quantitatively. 
The  change  in  flame  liftoff  height  shown  in  Tables  4  and  5,  and 
Fig.  11  indicate  that  approximately  ten  times  the  amount  of 
02(a1Ag)  (approximately  5500  ppm)  was  needed  to  achieve  the 
same  enhancement  as  03  (approximately  500  ppm).  Since  the 
enhancement  of  flame  propagation  speed  was  established  with 
the  addition  of  03  (shown  in  the  part  I  companion  paper)  [44],  it 
can  be  used  to  quantify  the  enhancement  by  02(a1Ag).  The  C2H4 
laminar  and  lifted  flame  speed  enhancement  was  computed  using 
the  C2H4  kinetic  mechanisms  [48,49]  with  the  addition  of  the  03 
reactions  shown  in  part  I  of  this  work  [44].  The  laminar  and  lifted 
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Fig.  12.  Plot  of  computational  results  of  lifted  flame  speed  enhancement  with 
02(a1Ag)  and  03.  The  horizontal  error  bars  denote  the  propagation  of  uncertainty 
from  the  concentrations  found  experimentally,  while  the  vertical  error  bars  denote 
the  range  of  lifted  flame  speed  enhancement  when  using  the  different  published 
reaction  rates. 
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flame  speed  are  related  through  the  square  root  of  the  density  ratio 
[44],  therefore  the  percent  enhancement  of  the  flame  speeds  were 
comparable.  Furthermore,  all  of  the  experiments  were  performed 
for  large  flame  liftoff  heights  where  the  mixture  fraction  gradient 
is  small  and  there  is  the  closest  agreement  to  the  numerical  simu¬ 
lations  at  a  zero  mixture  fraction  gradient.  Therefore,  the  lifted 
flame  speed  enhancement  was  used  for  comparison  between 
experiments  and  numerical  simulations.  The  results  for  the  condi¬ 
tions  of  500  ppm  of  03  addition  at  3.61  kPa  and  6.73  kPa  showed 
approximately  1%  enhancement  of  the  lifted  flame  speed 
(Fig.  12).  Therefore,  since  the  flame  liftoff  height  change  for 
500  ppm  03  is  equivalent  to  5500  ppm  02(a1  Ag),  it  was  reasonable 
to  assume  that  their  effects  on  flame  speed  enhancement  were 
comparable.  The  overall  lifted  flame  speed  enhancement  would 
be  higher  because  of  the  kinetic-induced  hydrodynamic  enhance¬ 
ment  described  in  part  I  of  this  work  [44].  If  the  effect  of  03  at 
low  pressure  is  comparable  to  what  was  found  at  high  pressure, 
then  5500  ppm  of  02(a1  Ag)  will  give  approximately  2-3%  enhance¬ 
ment  of  the  lifted  flame  speed. 

3.3.  Computational  analysis  with  current  kinetics  of  02(a1Ag) 

Numerical  simulations  were  performed  in  order  to  explain  the 
enhancement  mechanism  with  the  addition  of  02(a1Ag).  The 
C2H4  combustion  mechanisms  [48,49]  with  the  03  reactions  added 
[44]  were  used  along  with  the  inclusion  of  02(a1  Ag)  reactions.  The 
reactive  and  collisional  quenching  rates  of  02(a1Ag)  with  oxygen 
and  inert  species  as  a  function  of  temperature  have  been  well  stud¬ 
ied  and  were  used  to  compute  the  cold  flow  transport  of  02(a1Ag) 
in  the  ICOS  system  and  prior  to  the  flame  in  the  combustion  sys¬ 
tem,  as  shown  in  Section  3.1.  These  02(a1Ag)  reactions  were  added 
to  the  C2H4/03  mechanisms.  The  reactions  of  02(a1  Ag)  with  the  fuel 
have  also  been  well  studied  at  298  K,  but  there  is  little  data  at 
intermediate  and  high  temperatures.  The  most  studied  02(a1Ag) 
reactions  applicable  to  a  combustion  system  are  for  H2-02  mix¬ 
tures.  These  reactions  of  hydrogen  and  oxygen  containing  species 
with  02(a1  Ag)  are  shown  in  Table  6  and  have  been  compiled  specif¬ 
ically  for  H2-02  mixtures  activated  by  plasma.  The  rates  were  ta¬ 
ken  from  previous  numerical  work  on  plasma-assisted  H2-02 
combustion  systems,  as  well  as  from  specific  reaction  rate  studies 
[27,29,50-53].  There  are  a  few  differences  that  arose,  specifically 
regarding  the  two  reactions 


Fig.  11.  Plot  of  experimental  results  of  flame  liftoff  change  with  O^a^g)  and  03 
concentration  for  a  plasma  power  of  80  W.  The  error  bars  denote  the  experimental 
uncertainties  in  the  concentrations  and  change  in  flame  liftoff  height. 


H  +  02  (a’Ag)  -  O  +  OH 
H2  +  02  (a1  Ag)  — >  H2+02 


(5) 

(6) 


T.  Ombrello  et  al.  / Combustion  and  Flame  157  (2010)  1916-1928 


1925 


Table  6 

Reaction  rates  of  02(a1Ag)  with  hydrogen  containing  species. 


Reaction 

Reaction  constant  (cm3/molecule/s) 

Temperature  dependence 

Activation  energy  (kj/mole) 

Reactive  quenching  of  02(al  A g) 

H  +  02  ->  O  +  OH  [50] 

5.00  x  10“9 

0 

60.3 

H  +  02(a1Ag)  ->  O  +  OH  [51] 

1.82  x  10“10 

0 

26.5 

H  +  O^a’Ag)  ->  O  +  OH  [29] 

6.55  x  10”11 

0 

21.0 

OH  +  02(a1Ag)  ->  O  +  H02  [51  ] 

2.16  x  lO”11 

0 

141.4 

OH  +  O^a’Ag)  -  H  +  03  [51] 

7.31  x  10“17 

1.44 

226.4 

H2  +  02(a1Ag)  ^  OH  +  OH  [51] 

2.82  x  10“9 

0 

141.4 

H2  +  02(a1Ag)  ->  H02  +  H  [51] 

4.13  x  10“12 

0 

151.6 

H20  +  02(a1Ag)  ->  OH  +  H02  [27] 

9.03  x  10“8 

0.5 

209.7 

H20  +  02(a1Ag)  ->  O  +  h202  [27] 

2.05  x  10“12 

0.5 

283.5 

Collisional  quenching  of  02(a1  A g) 

H2  +  O^a’Ag)  ->  H2  +  02  [52] 

2.16  x  10“13 

0 

21.6 

H2  +  02(a3Ag)  ->  H2  +  02  [29] 

1.68  x  10“12 

0 

32.0 

H  +  02(a1Ag)  ->  H  +  02  [51] 

6.97  x  10“16 

0 

0 

OH  +  O^a’Ag)  -»  OH  +  02  [51  ] 

5.65  x  10“18 

0 

0 

H02  +  02(a1Ag)  ->  H02  +  02  [51  ] 

5.65  x  10“18 

0 

0 

H20  +  02(a1Ag)  ->  H20  +  02  [51  ] 

5.65  x  10“18 

0 

0 

H202  +  02(a1Ag)  -  H202  +  02  [51  ] 

5.65  x  10“18 

0 

0 

Other  reactions  with  02(a*Ag) 

H  +  H02  ->  H2  +  02(a1Ag)  [51] 

3.32  x  10“12 

0 

2.4 

OH  +  O  ->  H  +  02(a1Ag)  [51] 

9.63  x  10“12 

0 

51.8 

03  +  OH  ->  H02  +  02(a1Ag)  [51  ] 

7.97  x  10“13 

0 

8.3 

03  +  H02  ->  OH  +  02  +  02(a1Ag)  [51  ] 

1.66  x  10“14 

0 

8.3 

H02  +  H02  H202  +  02(a1Ag)  [51  ] 

1.49  x  lO”11 

0 

4.2 

H202  +  O  ->  H20  +  O^a’Ag)  [51  ] 

6.97  x  10“13 

0 

17.7 

There  was  a  small  difference  in  reaction  (5)  between  Refs.  [29]  and 
[51  ],  as  well  as  in  reaction  (6)  between  Refs.  [29]  and  [52],  as  shown 
in  Table  6.  The  differences  in  the  reported  reaction  rates  will  be  dis¬ 
cussed  in  relation  to  the  sensitivity  of  the  computed  flame  speed 
enhancement  shortly. 

The  reactions  were  added  to  the  C2H4/03  kinetic  mechanisms 
and  the  results  are  shown  in  Fig.  12  using  the  concentrations  of 
03  and  ([^(a1  Ag)  from  the  experiments.  The  results  of  enhancement 
using  the  02(a1Ag)  concentrations  found  in  the  experiments 
showed  flame  speed  enhancement  of  more  than  5%.  When  com¬ 
pared  to  the  results  of  lifted  flame  speed  enhancement  by  03 
(which  was  found  to  be  approximately  1%)  there  was  a  large  dis¬ 
crepancy  compared  with  the  experimental  results,  which  was  well 
outside  the  uncertainties  for  the  system.  For  example,  in  the  exper¬ 
iments  at  3.61  kPa,  the  change  in  flame  liftoff  height  for  500  ppm 
03  was  approximately  equal  to  5500  ppm  02(a1Ag)  (Fig.  11). 
According  to  the  numerical  simulations,  there  was  more  than  a  fac¬ 
tor  of  five  difference  in  the  enhancement  of  flame  speed.  Consider¬ 
ing  that  the  03  and  02(a1Ag)  should  both  enhance  the  lifted  flame 
speed  similarly  with  the  lean  and  rich  enhancement  more  than 
stoichiometric  [31,32],  there  is  a  significant  error  in  the  02(a1Ag) 
kinetic  calculations.  With  regard  to  the  differences  in  the  rates 
for  reactions  (5)  and  (6),  the  deviations  in  the  flame  speed 
enhancement  was  no  more  than  ±4%.  The  vertical  error  bars  in 
Fig.  12  show  the  negligible  enhancement  differences  when  using 
the  two  published  rates  for  reactions  (5)  and  (6)  and  therefore 
the  sensitivity  of  the  flame  speed  to  the  differences  in  reaction 
rates  do  not  account  for  the  significant  deviations  from  the 
experiments. 

To  understand  what  caused  the  significant  flame  speed 
enhancement  by  02(a1Ag),  a  rate  of  production  plot  of  02(a1Ag) 
in  the  flame  is  shown  in  Fig.  13.  The  major  consumption  pathways 
of  02(a1Ag)  was  from  the  branching  reaction  with  H,  with  some 
contribution  of  collisional  quenching  by  H2.  The  reaction  of 
02(a1Ag)  with  H  will  significantly  enhance  the  flame  since  it  is  a 
primary  radical  branching  reaction.  Therefore,  it  was  reasonable 
that  the  enhancement  was  so  large  in  the  numerical  simulations. 

The  changes  in  the  radical  pool  in  the  earlier  stages  of  the  flame 
were  plotted  and  are  shown  in  Fig.  14.  In  the  earlier  stages  of  the 
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Fig.  13.  Rate  of  production  plot  of  02(a1Ag)  superimposed  on  the  temperature 
profile  showing  the  major  consumption  pathways  of  02(a1Ag)  with  current 
published  rate  data  with  hydrogen  containing  species. 


flame  where  the  temperature  is  slightly  elevated  between  400  K 
and  500  K,  the  02(a1  Ag)  begins  to  be  consumed,  causing  a  decrease 
in  the  C2H4  concentration  and  a  subsequent  increase  in  O  and  OH. 
The  increase  in  concentration  of  OH  would  provide  chemical  heat 
release  through  subsequent  reactions  earlier  in  the  flame  to  en¬ 
hance  the  overall  flame  speed  significantly.  These  results  of  higher 
concentrations  of  OH  in  the  earlier  stages  of  the  flame  leading  to 
chemical  heat  release  and  enhanced  flame  speed  was  shown  in 
part  I  of  this  work  through  the  study  of  03  addition  [44].  The  in¬ 
crease  in  the  radical  pool  concentration  of  O  from  the  reactions 
with  02(a1Ag)  were  investigated  in  more  depth  through  a  rate  of 
production  analysis.  The  results  showed  that  the  major  pathway 
for  O  consumption  was  from  the  reaction  with  the  parent  fuel, 
C2H4,  and  its  fragment,  CH3.  The  end  results  of  flame  propagation 
enhancement  came  from  the  increased  radical  pool  concentration 
and  extraction  of  chemical  heat  release  earlier  in  the  flame  com¬ 
pared  to  the  results  with  no  02(a1Ag)  addition.  Fig.  15  clearly  indi- 
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Fig.  14.  Temperature  and  mole  fraction  profiles,  dashed  lines  =  w/o  O^a^g),  solid 
lines  =  w/02(a1Ag). 


Fig.  15.  Heat  release  versus  temperature  showing  elevated  levels  of  chemical  heat 
release  at  lower  temperatures  with  02(a1Ag)  addition. 

cates  the  elevated  levels  of  chemical  heat  release  result  by  showing 
the  total  volumetric  heat  release  as  a  function  of  temperature  in 
the  flame.  There  is  significantly  more  heat  release  between  800  K 
and  1500  K  with  02(a1Ag)  addition,  which  aligns  well  with  the 
peak  consumption  of  02(a1Ag)  and  elevated  radical  concentrations 
shown  in  Figs.  13  and  14,  respectively. 


3.4.  Inclusion  of02(a7Ag)  quenching  by  hydrocarbons 

With  an  understanding  of  how  02(a1Ag)  enhanced  the  flame 
speed  in  the  numerical  simulations,  it  becomes  apparent  that  there 
are  two  possible  explanations  for  the  discrepancy  shown  in  Fig.  12 
with  respect  to  the  enhancement  by  03  and  02(a1  Ag).  First,  a  signif¬ 
icant  concentration  of  the  02(a1Ag)  could  quench  before  reacting 
with  H.  The  collisional  quenching  reactions  of  02(a1Ag)  with  the 
parent  fuel  or  its  fuel  fragments  could  decrease  the  concentration 
significantly,  therefore  leading  to  less  enhancement.  To  achieve 
the  same  enhancement  wih  02(a1Ag)  as  the  calculations  with  03, 
there  would  need  to  be  approximately  900  ppm  02(a1Ag)  (shown 
in  Fig.  12  for  3.61  kPa).  Therefore,  approximately  4600  ppm  of 
02(a1Ag)  would  have  to  collisionally  quench,  allowing  only  a  small 
fraction  of  the  original  concentration  to  react  with  H.  The  enhance¬ 
ment  would  be  less,  especially  if  the  dominating  reactions  of 
02(a1Ag)  are  with  the  parent  fuel  and  not  providing  chain  branch¬ 
ing  as  in  the  reaction  of  H  with  02(a1Ag).  Second,  there  could  be  a 
combination  of  reactive  and  collisional  quenching  pathways  for 


Table  7 

Reaction  rates  of  02(a1Ag)  with  hydrocarbon  species.  (How)  =  estimated  rate  with 
low  activation  energy,  (*high)  =  estimated  rate  with  high  activation  energy,  (*exp. 
fit)  =  estimated  rate  with  activation  energy  to  fit  trend  of  experimental  results, 
(*)  =  estimated  rate. 


Reaction 

Reaction 

Temperature  Activation 

constant  (cm3/  dependence 

energy 

molecule/s) 

(1<J) 

CH4  +  O^a’Ag)  ->  CH4  +  02  [42] 

1.40  x  10-18 

0 

0 

CH4  +  O^a’Ag)  -►  CH3  +  H02  [54] 

6.14  x  10“12 

0 

149.0 

C2H4  +  02(a1Ag)  ->  C2H4  +  02  (How) 

7.71  x  10“16 

0 

15.0 

C2H4  +  02(a1Ag)  ->  C2H4  +  02  (*high)  3.12  x  10“13 

0 

30.0 

C2H4  +  02(a1Ag)  -»  C2H4  +  02 

5.46  x  10-10 

0 

48.6 

(*exp.  fit) 

C2H4  +  OzCa^g)  ->  C2H3  +  H02  (*) 

7.01  x  10  11 

0 

146.5 

02(a1Ag)  that  could  be  responsible  for  the  enhancement  observed 
in  the  experiments. 

The  reactive  and  collisional  quenching  rates  for  some  of  the 
hydrocarbon  species  were  added  to  the  kinetic  mechanism  one  at 
a  time  to  test  the  sensitivity  of  flame  speed  enhancement.  In  Ta¬ 
ble  7  is  a  list  of  the  reactions  along  with  their  rates.  Initially  the 
reactions  with  CH4  of 

CH4  +  02  (a1  Ag)  — >  CH4  +  02  (7) 

CH4  +  02  (a1  Ag)  -►  CH3  +  HOz  (8) 

were  added,  but  did  not  change  the  flame  speed  enhancement, 
which  was  reasonable  considering  the  low  concentrations  of  CH4 
in  the  system.  Next,  noting  that  the  inclusion  of  02(a1  Ag)  collisional 
quenching  by  the  parent  fuel  H2  in  an  H2-02  system  via  reaction  (6) 
was  found  in  previous  numerical  investigations  to  be  significant 
and  decreased  the  effectiveness  of  02(a1Ag)  enhancement  [29],  it 
is  reasonable  to  include  collisional  quenching  by  the  parent  fuel  via 

C2H4  +  02(a1Ag)  C2H4  +  02  (9) 

for  our  hydrocarbon  fueled  combustion  system.  There  are  some  col¬ 
lisional  quenching  rates  for  02(a1Ag)  with  C2H4  and  other  small 
hydrocarbon  fuels,  but  they  are  only  at  298  K  [42,52].  To  the  authors 
knowledge,  there  are  no  verified  quenching  rates  of  hydrocarbon 
species  (specifically  C2H4)  with  02(a1  Ag)  in  the  intermediate  to  high 
temperature  range  where  they  would  be  the  most  important  for  a 
flame  system. 

The  work  of  Borrell  and  Richards  found  that  the  temperature 
dependence  of  02(a1Ag)  quenching  by  H2  was  approximately 
Arrhenius  [53],  and  that  other  species  also  follow  an  Arrhenius 


Fig.  16.  Arrhenius  temperature  dependent  collisional  quenching  rates  of  02(a1Ag) 
by  H2  and  C2H4.  The  rates  with  C2H4  are  estimates. 
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temperature  dependence.  Therefore,  it  was  assumed  that  C2H4 
might  also  follow  this  temperature  dependence  for  quenching 
02(a1Ag).  Originating  with  the  quenching  rate  of  02(a1Ag)  by 
C2H4  at  298  K  [42,52],  an  Arrhenius  temperature  dependent  rate 
was  estimated  to  explain  the  trend  discrepancy  shown  in  the 
experimental  results  of  Fig.  12.  Three  different  temperature  depen¬ 
dencies  were  chosen  with  a  range  of  activation  energies.  Previously 
published  temperature  dependent  collisional  quenching  rates  of 
02(a1Ag)  have  shown  that  the  activation  energy  range  is  typically 
between  15  kj/mole  and  20  kj/mole  [55].  Furthermore,  the  activa¬ 
tion  energy  for  the  collisional  quenching  of  02(a1Ag)  by  H2  is  as 
high  as  32  kj/mole  [29].  Therefore,  the  range  of  15-30  kj/mole 
was  chosen  for  the  activation  energies  of  reaction  (9),  with  a  reac¬ 
tion  constant  chosen  in  order  to  agree  with  published  rates  at 
298  K  [42,52].  The  envelope  of  reaction  rates  for  (9)  using  activa¬ 
tion  energies  from  151<J/mole  to  30  kj/mole  is  shown  in  Fig.  16 
and  is  labeled  as  “Quenching  A”. 

Computations  were  performed  using  the  rates  within  the  enve¬ 
lope  of  “Quenching  A”  in  Fig.  16  and  the  results  are  shown  in 
Fig.  17.  The  flame  speed  enhancement  decreased  slightly,  but  not 
enough  to  explain  the  discrepancy.  In  an  attempt  to  explain  the 
disagreement,  a  rate  for  reaction  (9)  was  chosen  in  order  to  bring 
the  enhancement  by  02(a1Ag)  down  to  the  enhancement  by  03. 
To  accomplish  this,  a  reaction  constant  of  5.46  x  10-10  cm3/mole- 
cule/s  and  an  activation  energy  of  48.6  kj/mole  were  chosen  with 
the  temperature  dependence  shown  in  Fig.  16  as  “Quenching  B.” 
By  using  the  “Quenching  B”  rate  for  reaction  (9),  the  results  of 
enhancement  by  03  and  02(a1Ag)  were  approximately  equal, 
therefore  agreeing  with  the  trends  of  the  experimental  results 
(Fig.  17).  A  rate  of  production  analysis  of  02(a1Ag)  was  performed 
with  the  “Quenching  B”  rate  and  the  major  consumption  pathway 
shifted  to  reaction  (9)  with  negligible  consumption  by  reaction  (5). 
Additionally,  the  concentration  profile  of  02(a1Ag)  showed  a  more 
rapid  decrease  in  the  earlier  stages  of  the  flame,  but  no  appreciable 
increase  in  C2H4  consumption  or  O  and  OH  production.  Further¬ 
more,  the  quenching  pathway  involves  electronic-to-vibrational- 
translational  energy  transfer  [56]  which  releases  so  little  energy 
that  the  translational  temperature  changed  negligibly.  The  com¬ 
puted  temperature  profiles  confirmed  the  negligible  increase  in 
temperature,  and  therefore  that  a  significant  concentration  of 
02(a1  Ag)  was  consumed  and  did  not  affect  the  flame  in  the  process. 
The  high  activation  energy  and  hence  strong  temperature  depen¬ 
dence  of  reaction  (9)  given  by  the  estimated  rate  “Quenching  B” 
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Fig.  17.  Plot  of  computational  results  of  lifted  flame  speed  enhancement  with 
02(a1Ag)  and  03  using  the  estimated  collisional  quenching  rate  of  C2H4  with 
02(a1Ag)  from  Fig.  14.  “Quenching  A”  =  inclusion  of  temperature  dependent 
quenching  of  02(a1Ag)  by  C2H4  with  Ea  =  30  kj/mol,  “Quenching  B”  =  inclusion  of 
temperature  dependent  quenching  of  02(a1Ag)  by  C2H4  to  fit  trend  of  experimental 
results.  The  error  bars  denote  the  variation  in  lifted  flame  speed  enhancement  when 
using  different  C2H4  kinetic  mechanisms. 
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mitigated  the  computed  enhancement  discrepancy,  but  the  rate 
appears  unreasonably  high  and  has  not  been  validated  and  there¬ 
fore  requires  further  investigations. 

Therefore,  the  last  possible  explanation  of  the  trend  discrepancy 
in  computed  enhancement  lies  in  the  reactive  quenching  of 
02(a1Ag)  by  C2H4  and  its  fragments.  The  first  assumption  was  to 
decrease  the  activation  energy  of  the  reaction 

C2H4  +  02(a1Ag)  ->  C2H3+H02  (10) 

by  the  energy  contained  within  02(a1  Ag).  This  equates  to  decreasing 
the  activation  energy  by  0.98  eV  (94.5  kj/mole),  and  the  rate  is 
shown  in  Table  7.  The  inclusion  of  this  reaction  in  the  kinetic  mech¬ 
anism  did  not  result  in  any  change  in  flame  speed  enhancement  be¬ 
cause  the  rate  is  slow  in  comparison  to  other  reaction  pathways 
with  C2H4  and  02(a1Ag).  Beyond  reaction  (10)  there  could  be  other 
possible  product  pathways  which  have  been  investigated  through 
quantum  calculations  with  C2H4  [57,58].  The  pathways  show  that 
02(a1Ag)  can  attack  the  double  carbon  bond  to  split  the  parent  fuel 
molecule,  possibly  providing  significant  enhancement  of  fuel  oxida¬ 
tion  rates  by  producing  CH20  and  other  hydrocarbon  fragments,  but 
the  rates  are  not  known. 

Therefore,  the  results  suggest  that  the  probable  reactive  scheme 
has  both  collision  and  reactive  quenching  of  02(a1  Ag)  by  the  parent 
fuel  and  its  fragments.  There  remain  many  unknowns  as  to  the  ki¬ 
netic  mechanisms  for  02(a1Ag)  with  hydrocarbons  under  flame 
conditions.  The  lack  of  rate  data  for  these  reactions  in  the  interme¬ 
diate  to  high  temperature  range  which  is  applicable  to  combustion 
systems  remains  a  significant  obstacle  and  requires  further  inves¬ 
tigation.  Nevertheless,  the  results  from  this  investigation  have  pro¬ 
vided  the  first  experimental  data  set  of  flame  propagation 
enhancement  by  02(a1Ag)  which  provides  a  foundation  for  future 
investigations. 


4.  Conclusions 

The  present  study  provides  a  promising  approach  to  isolate  and 
transport  plasma-produced  excited  species  for  the  kinetic  study  of 
plasma-assisted  combustion.  By  separating  the  plasma  and  com¬ 
bustion  system,  specific  plasma-produced  species  can  be  isolated 
and  measured  while  minimizing  the  complications  of  other  plas¬ 
ma-flame  interactions.  Isolation  of  the  specific  effects  of  individual 
plasma-produced  species  will  have  a  significant  impact  on  the 
development  of  detailed  plasma-flame  kinetic  mechanisms. 
Through  the  current  work  a  platform  to  study  quantitatively  the 
enhancement  effects  of  plasma-produced  02(a1Ag)  on  C2H4  lifted 
flame  propagation  speeds  was  developed.  It  was  found  quantita¬ 
tively,  for  the  first  time,  that  02(a1  Ag)  enhances  flame  propagation. 
The  addition  of  NO  to  the  plasma  afterglow  allowed  for  an  order  of 
magnitude  increase  in  the  02(a1Ag)  concentration  at  a  given  resi¬ 
dence  time  by  removing  the  quenching  species  of  03  and  O.  The 
NO  was  extremely  effective  because  of  the  catalytic  cycle  to  re¬ 
move  03  and  O,  as  well  as  having  a  negligible  effect  on  flame  speed. 
The  02(a1Ag)  was  produced  in  concentrations  of  over  5000  ppm  to 
enhance  flame  propagation  of  C2H4  lifted  flames  by  several  percent 
at  3.61  kPa  and  6.73  kPa. 

Numerical  simulations  using  currently  published  collisional  and 
reactive  quenching  reactions  have  shown  that  there  is  a  significant 
discrepancy  in  the  predicted  enhancement  compared  with  the 
trends  found  in  the  experiments.  The  pathways  of  enhancement 
found  in  the  simulations  showed  that  the  branching  reaction  of 
02(a1Ag)  with  H  provided  O  and  OH  early  in  the  reaction  zone 
and  increased  chemical  heat  release  and  flame  propagation 
enhancement.  The  lack  of  temperature  dependent  quenching  rate 
data  of  02(a1Ag)  by  hydrocarbon  species  was  suspected  to  be  the 
main  cause  for  the  discrepancy.  Estimations  of  the  temperature 
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dependent  collisional  quenching  rate  of  02(a1Ag)  by  C2H4  has 
shown  good  agreement  with  the  experimental  results  trends,  but 
the  suggested  rate  appears  to  be  unreasonably  high  and  needs  to 
be  validated  in  future  studies.  Furthermore,  the  reactive  quenching 
pathways  and  their  rates  remain  unknown.  Therefore,  a  combina¬ 
tion  of  both  the  collisional  and  reactive  quenching  rates  of 
C^^Ag)  with  hydrocarbon  species,  specifically  the  parent  fuel 
molecule,  are  required  in  order  to  correctly  explain  combustion 
enhancement.  The  experimental  results  have  therefore  provided 
the  first  data  of  the  isolated  effect  of  02(a1Ag)  under  flame  condi¬ 
tions,  which  is  of  paramount  importance  for  the  development  of 
reaction  pathways  and  plasma-flame  kinetic  mechanisms.  Future 
investigations  are  being  targeted  at  establishing  temperature 
dependent  quenching  rates  of  O^a^g)  to  enable  more  accurate 
and  predictive  modeling  of  the  plasma-flame  interaction. 

The  experimental  results  also  imply  that  when  energy  is  cou¬ 
pled  into  specific  plasma-produced  species,  there  is  no  require¬ 
ment  for  the  control  of  heat  loss.  If  energy  is  coupled  into  a 
reactive  flow  to  raise  the  translational  gas  temperature  only,  there 
needs  to  be  careful  thermal  management.  Whereas,  when  energy  is 
coupled  into  energy  modes  of  specific  species,  no  thermal  manage¬ 
ment  needs  to  be  considered.  The  energy  contained  within  the  spe¬ 
cies  can  be  transported  for  significant  distances  and  residence 
times  and  extracted  at  the  combustion  reaction  zone.  Of  interest 
to  note  is  that  the  plasma  power  remained  constant  at  80  W  for 
all  conditions  with  03  and/or  02(a1Ag)  in  the  current  experiments. 
Therefore,  the  energy  coupled  into  the  flow  by  the  plasma  was 
recovered  at  the  flame.  More  02(a1Ag)  was  produced  than  03  for 
the  same  plasma  power  because  of  the  difference  in  energy  to  pro¬ 
duce  these  species. 

Lastly,  in  and  immediately  downstream  of  an  oxygen  plasma 
there  would  be  both  O  and  02(a1Ag),  but  no  03  since  it  is  the  prod¬ 
uct  of  a  time-dependent  recombination  of  O  and  02.  The  0  atoms 
will  most  likely  be  more  reactive,  but  there  will  be  higher  concen¬ 
trations  of  02(a1Ag).  Therefore,  if  oxygen  plasma  is  located  closer 
to  a  combustion  reaction  zone,  the  effects  of  02(a1Ag)  can  become 
more  pronounced. 
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